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Abstract. We introduce a long wave scaling for the Vlasov-Poisson 
equation and derive, in the cold ions limit, the Korteweg-De Vries equa- 
tion (in ID) and the Zakharov-Kuznetsov equation (in higher dimen- 
^ | sions, in the presence of an external magnetic field). The proofs are 

D , based on the relative entropy method. 



1. The long wave scaling of the Vlasov-Poisson equation 

1.1. The Vlasov-Poisson system for ions with small mass electrons. 

We consider the Vlasov-Poisson system which describes the evolution of ions 
in a plasma. We assume that due to their small mass, electrons in the plasma 
instantaneously reach their thermodynamic equilibrium, so that their density 
n e follows the Maxwell-Boltzmann law: 

(1.1) n e = e*/ T °, 

denoting by (ft the electric potential and T e > the scaled temperature of 
electrons, taken equal to 1 in the following. We then obtain the Vlasov- 
Poisson system for ions (where t > 0, x 6 T d or R d , v G M. d , d = 1, 2, 3): 

' dtf + v • V x f + (E + v A b) ■ V v f = 0, 
E = -V x( ft, 

(1.2) 



-A x (ft + e*= / fdv, 
_ f\t=o = fo- 

We refer to [26] for a more complete discussion on this system. In these equa- 
tions, f(t,x,v) stands for the distribution function of the ions, which allows 
to describe their statistical distribution: this means that f(t,x,v)dxdv ex- 
presses the density of ions, at time t, whose position is close to x and whose 
velocity is close to v. Furthermore, with usual notations, E is the electric 
field (generated by electrons and ions themselves) , while b is a fixed external 
magnetic field. 

It is interesting to notice that for a very particular class of (singular) data, 
namely Dirac functions in velocity, which we shall call monokinetic data, 

(1-3) f(t,x,v) = p{t,x)5 v=u{ t )X) , 

where 5 denotes the Dirac delta function, we get that / is a solution in the 
sense of distributions to (|1.2p (with some relevant electric field) if and only 
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if (p, u) satisfies the following hydro dynamic system, which corresponds to 
the pressureless Euler-Poisson system: 

' dtp + V x ■ (pu) = 0, 
dtu + u ■ V x u = E + u A b, u = (ui, u 2 , u 3 ), 
E = -V x </>, 
- A x (f) + = p. 



(1.4) 



From the physical point of view, this corresponds to the cold ions assump- 
tion, that corresponds to the situation where ions have zero (kinetic) tem- 
perature: any function of the form (jl.3p indeed satisfies 



(1.5) 



f(t, x, v)\v — u(t, x)\ 2 dv = 0. 



We also mention that a standard model in plasma physics is obtained after 
linearizing the Maxwell-Boltzmann law, which yields from (jl.2p the following 
system: 

f dtf + v ■ V x f + E ■ V v f = 0, 
E = -V T A 



;i.6) 



■ + 



fdv-1, 



k f\t=o — fo- 

This linearization produces some considerable simplifications in the analysis. 

1.2. The long wave scaling. In some recent independent works, Lannes, 
Linares and Saut |29) in the first hand, and Guo and Pu \22\ [38] in the 
other hand, have rigorously studied the long wave limit of the pressureless 
Euler-Poisson system. Precisely, one looks for solutions to (|l,4p under the 
form: 

' p = 1 + epi{£ 1/2 {xi - t),e 1/2 x 2 ,e 1/2 x 3 , e 3/2 t) + e 2 p 2 + 
<f> = e0i(e 1/2 (x! - t),e 1/2 x 2 ,e 1 / 2 x 3 , e z/2 t) + e 2 4> 2 + 

,( 2 ) 



(1.7) 



Ul = e 1 / 2 u^(e 1 / 2 (x 1 -t),e 1 / 2 x 2 ,e 1 / 2 x 3 ,e^ 2 t) + eu {2 ' + 



u 2 = e 1 / 2 u i 2 1 \e 1 / 2 (x 1 -t),e^' z x 2 ,e^' z x 3 ,e^' z t) + eu^> + 



,1/2 



1/2, 



r 3/2. 



,( 2 ) 



( u 3 = e 1 / 2 u < i\e l / 2 (x 1 -t),s 1 / 2 x 2 ,s 1 / 2 x 3 ,e 3 / 2 t) + ou , 

This leads to the study of the behaviour, as the parameter e goes to 0, of 
the solutions to the rescaled system: 

( ed t p - 8 Xl p + V x ■ ((1 + ep)u) = 0, 

edtu — d xi u + eu ■ V x u = E + e^ 1 ^ 2 u A b, u = (ui,u 2 ,u 3 ), 
E = -V x (j), 

- e 2 A x cf> + = 1 + ep, 



(1.8) 



In the limit, one expects to obtain some (non-linear) dispersive equations, as 
explained in the papers by Zakharov-Kuznetsov |45j and Laedke-Spatschek 



FROM VLASOV-POISSON TO KDV AND ZK 



3 



[28] . Indeed, in ID, both [29] and |22j establish the convergence to the 
Korteweg-de Vries (KdV) equation, which we recall below: 

(1.9) ^ + 01^01 + 5^01 = 0. 

In 2D and 3D, which corresponds to the setting of (|1.8p . [29j and [38j give 
the derivation of a higher dimensional generalization of the KdV equation, 
which is referred to as the Zakharov-Kuznetsov (in short ZK) equation: 

(1.10) a t <j> 1 + t f> 1 a x <i>x + d xl Au 1 = o. 

With another kind of anisotropic (in space) scaling, Pu has also derived 
in 2D the Kadomstev-Petviashvili II (in short KP-II) equation |38| : 

(1.11) d Xl (d t <f>i + (pid Xl (pi + dl lXlXl 4>i) + dl 2X2 cpi = 0. 

Over the past years, there have been many mathematical studies of long 
wave limits towards KdV (and higher dimensional generalizations). Let us 
cite some works (this list is not meant to be exhaustive) concerning the 
following PDEs: 

• Nonlinear Schrddinger (including Gross-Pitaevskii) equations: Chi- 
ron and Rousset |14j . and Bethuel, Gravejat, Saut and Smets [U [7] 
(with different methods), 

• General Hyperbolic systems: Ben-Youssef and Colin [3], Ben-Youssef 
and Lannes [5], 

• Water Waves: Craig [16], Schneider and Wayne |43J, Bona, Colin and 
Lannes [8], Alvarez-Samaniego and Lannes [2] , Duchene |17| 118] . 

and many others. 

The fact that there exist very singular solutions to the Vlasov-Poisson 
system (|1.2p which precisely yield the pressureless Euler-Poisson system (ll.4p 
suggests that is should also be possible to study (|1.2p in a long wave regime. 
Following this idea, we would like to look for solutions of the form: 



1.12) 



f £ (t, x, v) = e d f £ (s 3 ' 2 t, e 1 ' 2 (x 1 - t),e 1 / 2 x 2 , e^ 2 x 3 , e^v), 
<f> e {x, v) = e0 £ (e 3/2 t, e 1/2 (xi - t),e 1 / 2 x 2 ,e l/2 x z ). 



The normalization is chosen in order that the scaling of the two first moments 

f f v d 

p £ := f f £ dv and u £ := jj dv matches with the Ansatz in (|1.7p . 

At some point, we will also have to somehow impose that the function f £ 
is "close" to a Dirac function, in order to reach the Euler-Poisson dynamics. 

Therefore, we propose a long wave scaling for the Vlasov-Poisson equation 
(in dimension d = 3), which we introduce now: 



(1.13) 



i=e'l 2 t, x = e 1 / 2 x, v = e x v, 
f(t,x,v) = e~ 3 f{t,x,v), 
4>(t, x) = £~ l (j)(t, x), 
E{i,x) =e- 3/2 E(t,x). 
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With a slight abuse of notation (we forget the tildes), we obtain the rescaled 
equations: 

' E^ 2 d t j + e 3 / 2 v ■ V x f + e- 1 (e^E + ev A ft) • V v f = 0, 

E = 



(1.14) 



e 2 A x + e E 



fdv, 



. /|t=o — /o- 

Finally, there only remains to perform the shift with respect to the first 
spatial variable: 

' x\ = x\ — t, 

f(t,x 1 ,x 2 ,x 3 ,v) = f(t,x,v), 

4>(t,x 1 ,x 2 ,x 3 ) = (j)(t,Xl,X2,X 3 ), 

^ E(t,x 1 ,x 2 ,x 3 ) = E(t,X!,x 2 ,x 3 ). 

With another abuse of notation, we end up with the rescaled Vlasov-Poisson 
system: 



(1.15) 



(1.16) 



ed t f £ -d Xl f e + evV x f e +[E e + 



v A b 



v„/ 6 = o, 



-V, 



e 2 A x (f> £ + e £ 



„ /e,|t=0 — /e,0- 

To the best of our knowledge, although this scaling seems very natural, 
it has never been introduced yet, even in the physics literature. Our goal 
in this work is to study the behaviour, as e — > 0, of solutions to (|1.16p . 
that also get close (in some sense that we shall precise later) at initial time 
to monokinetic data. According to the previous discussion, it is natural to 
expect to obtain the ZK equation in the limit. 

The relations between the different systems are summarized in the follow- 
ing diagram: 



Vlasov- Poisson 


Cold ions limit 


Pressureless Euler-Poisson 





Combined 
cold ions and long wave limit 



Long wave limit 



KdV or ZK 



Let us briefly comment on this picture. As already explained, the link be- 
tween Vlasov-Poisson and Pressureless Euler-Poisson is given by monokinetic 
type of data. Unfortunately, these data are way too singular to fit in the 
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known Cauchy theories (although some results for measure data are actually 
available in ID, we refer to the work of Zheng and Majda |46j). As a matter 
of fact, we are not even aware of any result proving rigorously the conver- 
gence to the pressureless Euler-Poisson system for data that would converge 
in some sense to monokinetic data; the stability estimates that would be 
needed are indeed missing. 

As already said, the long wave limit from the Pressureless Euler-Poisson 
system has been performed in |29 | 122 } [38] . One important step is to build a 
solution in an interval of time which is independent of s. In order to study 
the long wave limit of the Vlasov-Poisson equation, instead of trying to derive 
the Pressureless Euler-Poisson system, the idea is to perform simultaneously 
the cold ions and long wave limits. We shall start from global weak solutions 
to the Vlasov-Poisson equation, and therefore we will not have to face the 
difficulty of finding uniform lifespans. 

To prove such a result, we shall rely on a classical energy method, namely 
the relative entropy method. The idea originates in the work of Yau |44| 
on the hydrodynamic limit of some Ginzburg-Landau equation. It was inde- 
pendently brought in kinetic theory by Golse in |10| (in the context of the 
incompressible Euler limit of the Boltzmann equation) and by Brenier in |13) 
(in the context of the quasineutral limit of the Vlasov-Poisson equation). 

The basic principle of the relative entropy strategy consists in modulating 
some well-chosen functional that has to be conserved or dissipated by the 
physical system (for instance, the good choice is the entropy for the Boltz- 
mann equation). The modulation is obtained in terms of the solution to 
the target equation. One has to ensure that this new functional allows to 
"measure" in a certain sense the distance between the solution to the original 
system and that to the target equation. Then one has to prove that the func- 
tional that has been constructed is a Lyapunov function for the system: this 
follows from exact computations and algebraic identities. The computations 
can be more or less lengthy and tedious. It can be also very technical to 
justify these. 

The estimates which can be obtained have to be understood as stability 
estimates: for instance, the results proved in this work can be interpreted as 
the stability of monokinetic data, in the long wave regime, with a dynamics 
dictated by some the KdV or ZK equations. This also strongly suggests that 
for our long wave limit, there are stability phenomena (a la Penrose) at stake, 
exactly like for the case of the quasineutral limit (the effects of instabilities 
for the latter limit are briefly discussed for instance in |21|): more precisely, 
two stream instabilities are bound to destabilize the system and make the 
long wave limit fail (but these are avoided when one considers monokinetic 
data). On the topic of instabilities for the Vlasov-Poisson system, we refer 

to [211 E31 ESI [301 ED EZl E21 E31 • 

For the Vlasov equation, this method is remarkably well adapted to handle 
the cold ions limit, in other words the "convergence" to monokinetic data, 
as it will be clear later. This was observed for the study of the quasineutral 
limit of the Vlasov-Poisson system for electrons with fixed ions, as done 
by Brenier |13j (this was completed later by Masmoudi [36] , see also |20j). 



6 



DANIEL HAN-KWAN 



More recently, in |26j , we have studied the quasineutral limit of the Vlasov- 
Poisson for ions with small mass electrons (which corresponds to (ll.2p ). In 
that work, we have observed that this equation displays a LlogL structure 
that is reminiscent of that of the Boltzmann equation (we refer to the works 
of Saint- Raymond |40| I42| . |41] for the incompressible Euler limit, see also 
the recent paper of Allemand pQ); this will also play a crucial role in this 
work. 

It is worth noticing that the method provided in this paper can also be 
applied to study the KdV limit of the Euler-Poisson system, for data with 
only low uniform regularity. Indeed we can start from the global weak solu- 
tions built by Cordier and Peng |15j and use similar computations as in the 
present paper. 

To conclude this introduction, let us mention that Haragus, Nicholls and 
Sattinger in |27| relied on the KdV approximation of the Euler-Poisson sys- 
tem to study (formally and numerically) the interaction of solitary waves. 
It would be very interesting to start an analogous program for the Vlasov- 
Poisson equation. 

1.3. Organization of the paper. The paper is organized as follows: in 
Section [2j we provide the derivation of the KdV equation (see Theorem 12. 3p , 
starting from the ID Vlasov-Poisson equation with a linearized Maxwell- 
Boltzmann law. The exposure of this relatively simple case will allow us to 
lay down the basic principles of the relative entropy method applied to the 
long wave limit. In addition to this pedagogical interest, the existence of 
global solutions to the KdV equation allows to give stability estimates which 
are valid for all times. In Section [3l we present the main result of this paper, 
which is the derivation of the ZK equation, in Theorem 13.31 starting from 
the 3D Vlasov-Poisson equation with the full Maxwell-Boltzmann law. The 
proof will be much more technical, in particular due to the fact that only 
a LlogL type of control is available for the electric potential (instead of a 
L 2 bound which can be obtained with a linearized Maxwell-Boltzmann law) . 
We will need an unusually large number of correctors in the relative entropy. 
Finally, we give in two appendices some variants of our results (which can 
still be obtained with the relative entropy method) : in particular we present 
another scaling for the 2D Vlasov-Poisson system which yields the KP-II 
equation in the long wave limit. 

2. From the Vlasov-Poisson equation to the Korteweg-de 



In this section, we shall study the long wave limit of the ID Vlasov-Poisson 
system with a linearized Maxwell-Boltzmann law, that is (here (x, v) G T x 



Vries equation 



R) 



' e d t fe - d x f £ + e vd x f £ + E E d v f e = 
E e = -d x 4> £ , 



(2.1) 



< 




, fe,\t=0 — fefi- 
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2.1. Preliminaries. This system possesses an energy, which is conserved, 
at least formally: 

(2.2) S e (t) :=^J f E \v\ 2 dvdx + ~e J \d x <f> e \ 2 dx + ^ J \<P £ \ 2 dx. 

Using this energy, as well as the conservation of L XyV norms that can be 
obtained using the hamiltonian structure of the Vlasov equation, one can 
prove, following the work of Arsenev [3], the following theorem, which states 
the existence of global weak solutions to (|2.ip : 

Theorem 2.1. Let e > 0. Let f efi G L 1 n L°°(T x R) be a non-negative 
function such that the initial energy is bounded: 

(2.3) E e (0) :=^j f £fi \v\ 2 dvdx+^e J \d x <f> Bi0 \ 2 dx+~ J \^ \ 2 dx < +oo, 
where the initial electric potential <j) E Q is given by the elliptic equation: 

-e 2 d 2 x (f) Et0 + e(j) £fi = / / e o dv - 1. 



We also assume that: 

j fefidv dx = 1. 

Then there exists a non-negative global weak solution f e E L^°(L 1 Pi L°°(T x 
to (|2.ip . such that the energy is non-increasing: 



(2.4) Vt > t', £ £ (t) < S £ (t'), 

and such that the following local conservation laws for {p e := J f e dv,J £ : = 
J f £ vdv) are satisfied: 

(2.5) d t pe - ^dxPe + d x J £ = 0, 

(2.6) d t j £ - - £ d x j £ + d x (J \v\ 2 f e dv^j = -^d x (4> £ + i) 2 + e -d x \d x v £ \ 2 . 

For the KdV equation, that we recall below, 

(2.7) d t <h + <hd*<h. + d* xx <t>! = 0, 

we have in hand a famous global existence result which was proved in the 
seminal paper of Bourgain |llj : 

Theorem 2.2. The KdV equation is globally well-posed in H S (T) for s > 0. 
We will use this result only for large values of s. 

2.2. Formal derivation. It is quite enlightening to perform a formal analy- 
sis in this simple one-dimensional case, to understand how the KdV equation 
arises. Formally, one directly considers monokinetic data: 

(2-8) fe(t,X,v) = Pe(t,x)6 v=Ue ( t>x y 
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Then, as already explained, we obtain the following pressureless Euler-Poisson 
equation in a long wave scaling: 



(2.9) 



d t pe - ~d x p £ + d x (p £ u £ ) = 0, 



d t u £ - -d x u £ + u £ d x u £ 
e 

-dx<f>e, 



'-E £ , 



E 



P-: 



1, 



(2.10) 



We look for an approximate solution satisfying the Ansatz: 

p £ = 1 + epi + 0(e 2 ), 

</> e = fa+efo + 0(e 2 ), 

v u £ = u\ + eu 2 + 0(e 2 ). 

Plugging this Ansatz in the equations, and matching the different powers of 
e, we obtain a cascade of equations. 

• Conservation of charge equation: 

(2.11) 0(1): -d x p 1 + d x u 1 = 0. 

(2.12) O(e): d tPl - d x p 2 + d x { Pl ux) + d x u 2 = 0. 

• Momentum equation: 

(2.13) 0{e- v ) : -8^+8^ = 0, 

(2.14) 0(1): d t u 1 -d x u 2 + u l d x u 1 + d x ^ 2 = Q. 

• Poisson equation 

(2.15) 0{e): <f>i = pi, 

0{e 2 ): -d 2 xx 



(2.16) 



>i + (p2 



Pi- 



We clearly get from (|2TTj) and ([2TT5]) : 
(2.17) fa = pi = ui 

and in the other hand from (|2.14p and (|2.16p : 

<9 X (02 - u 2 ) = d x ((p2 - Pz) + d x (p 2 - u 2 ) 



(2.18) 



d xx 4>i + d t ui + 2uid x u\. 



Hence, from (|2.12p and (|2.18p . we conclude that <pi satisfies the KdV 
equation: 

(2.19) 2S^i + 301^01 + C^i = 0. 

2.3. Rigorous derivation. We now state the theorem which rigorously 
proves the convergence to KdV. 

Theorem 2.3. Let (/e,o)ee(o,i) ^ e a family of non-negative initial data sat- 
isfying the assumptions of Theorem \2.1\ and such that there exists C > 
with: 



(2.20) 



£ E (0) < C, Vee (0,1). 
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We denote by (/ £ ) e g(o,i) a family of non-negative global weak solutions to 
(|2,ip given by Theorem \2.1\ Let TL £ be the relative entropy defined by the 
functional: 
(2.21) 

Ue(t) := i J f e \v-ui-eu 2 \ 2 dvdx+^e J \d x fa - d x fa - ed x fa\ 2 dx 



where (m, fa, u 2 , fa) G [C([0, +oo[, i? s+2 (T))] 4 (with s > 3/2) satisfy the 
following system: 

( 2d t fa + 3fad x fa + d 3 xxx fa = 0, 

(2.22) I ui = fa, 

{ d x (u 2 - fa) = d t fa + fad x fa. 

Then there exist Ci,C 2 > 0, such that for all e G (0, 1), 

(2.23) Vt > 0, U £ (t) < H e (0) + I {CxUeis) + C 2 JI) ds. 

Jo 

Assuming in addition that there exists C 3 such that for any e G (0, 1): 

(2.24) n £ (0) < C 3V ^, 
then we obtain for any e G (0, 1): 

(2.25) Vt > 0, H £ {t) < C 3 e Clt VI + C ^' 1 ^. 

In order to get functions satisfying (|2.22|) . we can proceed as follows: 

• The existence of fa is ensured by Theorem 12.21 

• We accordingly set fa = u\. 

• The functions fa and u 2 can be seen as correctors. The last equation 
of (|2.22p allows to define them. Only the value of the function u 2 — fa 
is important. 

Remark 2.1. The estimate (|2.25p clearly implies that 

(2.26) \ [ f e \v - m - eu 2 \ 2 dv dx < C 3 e Clt ^/I + C 2 6 * ~ 1 y/i, 
2 J Ci 



and 

(2.27) ~ J (fa - fa - efaf dx < C 3 e Cit V^ + C- f^'^ yfl. 

From (|2.26p . denoting by f a weak limit of f £ , we deduce that necessarily 

(2.28) J f\v -m\ 2 dvdx = 0, 

which means that the limit temperature is zero (cold ions limit). 
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Remark 2.2. The estimate (|2,25p is valid for all times, but it is useful 
for times of order o(| loge|) (in other words, logarithmically growing times). 
This is slightly better than the times of validity obtained in the long wave 
limit of Euler-Poisson in [29 J, which are (translated in our framework) of 
order 0(1). 

From this theorem, we can deduce the following corollary: 

Corollary 2.1. Making the same assumptions as in the previous theorem, 
we obtain the weak convergences: 

' p £ ^ e ->o 1 in Lf'M 1 weak-*, 
(2.29) \ J e ^ e ^Q 0i in L^M 1 weak-*, 

k 4> e ^e^o 4>\ in LfM 1 weak-*. 
Proof of Corollary \ 2. 11 By conservation of the L 1 norm, we deduce that 

Ikllzi < C. 

Using this bound and the one coming from the energy inequality: 

f £ \v\ 2 dv dx < C, 

and by non- negativity of f £ , it is standard to deduce a L 1 control on J £ : 

\\PeWv < C. 

Therefore, there exist two non-negative measures p, J such that up to 
some extraction, we have: 

' Pe P, 

in the vague sense of measures. We have to show that p = 1 and J = (pi. 

We pass to the limit e — > (in the sense of distributions) in the Poisson 
equation: 

_ 2 . 

" '-'xx'+'e 

using the uniform L? bound on y/ed x <j) £ and <\> e (coming from the energy 
inequality), we deduce that we necessarily have p = 1, which proves the first 
claim. 

In the other hand, using the Cauchy-Schwarz inequality, we have 
(2.31) 

\J C - p e {ui +eu 2 )\ 2 {J fe(v - ux - eu 2 )dv) f 2 

< fe\V — Ul—EU2\ dv, 



(2.30) 



-£ 2 dl x (f) £ +e<f> E = p e -l, 



Pe f fedv 



|2 



The functional (p, J) — > f - — 6U2 " dx is convex and lower semi- 
continuous with respect to the weak convergence of measures (see [13J). As 
a consequence, the weak convergences in the vague sense of measures p e — 1 
and J t — J lead to: 

(2.32) f\J- Ul \ 2 dx < liminf f ^ ~ PeU ^ dx. 

J ^ i Pe 

By (|2.25p . we deduce that J = u\{= <pi). 
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To conclude, the uniqueness of the limit allows us to say that the weak 
convergences actually hold without any extraction. 

□ 

Remark 2.3. We can actually state strong convergence results. Indeed, 
in view of the preceeding proof, it is clear that (|2.26|) implies the "strong" 
convergence: 

2 



(2.33) 



J e - p e (m + eu 2 ) 



<C^ t V~e + C 2 — r 
L 2 (-1 



1 



VS. 



In the other hand, the control (|2.27l) means that <j> e converges strongly in 
L 2 to 4>i, as e goes to 0. 



2.4. Proof of Theorem l2.3l Relying on the fact the energy is non-increasing 
(and that T~L £ is built as a modulation of the energy), we have for all t > 
and all e G (0, 1): 

U&) = £ £ (t) + (H £ (t) - £ £ {t)) < £ £ (0) + (H e (t) - £ e (t)), 

which yields: 



u e (t) <n e (p) + J* J e t 



fe ( + zu 2 \ 2 - v(ui + eu 2 ) 



dv dx 



+ ej8 t 
+ fdt 



■\d x <j>i + ed x 4> 2 \ 2 - d x 4> £ {d x (t)i + ed x <p 2 ) 



dx 



1 



h + eh) 



dx ds 



:=n £ (0)+ l\h + I 2 + h)ds. 
Jo 



We are now going to study I\, I 2 and I3. The computations can be justified 
using only the local conservation laws fj2.5|) and (|2.6p . 



Study of I x . 

Using the fact that f £ satisfies the Vlasov-Poisson equation, we obtain the 
identity: 



J d t f £ ^\ui + eu 2 \ 2 - v(ui + eu 2 )^j dvdx 



f £ (ui + eu 2 - v) 



1 



1 



d x ui + vd x ui - d x u 2 + evd x u 2 



dv dx 



-p £ E £ u\ dx - \ p £ E £ u 2 dx. 



In order to obtain an hydrodynamic equation inside [...], we write: 



f £ (ui + eu 2 - v)(vd x ui) = J f E (ui + eu 2 - v){u\d x ui) 

-J fe\ui+eu 2 - v\ 2 d x U! + J ef E (u\ + eu 2 - v)u 2 d x ui. 
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After deriving with respect to time the term ^\u\ + eu2\ 2 — v(u\ + eu2), 
we get the following contribution in I±: 



d t ui + ed t u 2 



dv dx. 



f £ (ui + eu 2 - v) 

We now focus on the terms of order 0(l/e), for which we can write 
— J f £ (ui + EU 2 — v) — d x u\ dv dx 

= j f £ (ui+eu 2 -v)^(-d x ui + d x (pi)dv dx 



-p £ u\d x (f)i dx - j p £ u 2 d x (/) 1 dx + j -J £ d x <f>i dx, 



and 



- J ^p e E £ ui dx - J ^p £ u\d x <\)\ dx - J p £ uid x 4>2 dx 
In the other hand, we observe that we can write: 

/ 



x<P2 ~ Ux<, 



i dx. 



f £ (ui + eu 2 - v) 



d t ui + uid x ui - d x u 2 



dv dx 



I 



dv dx 



f £ (ui + eu 2 - v) d t ui + uid x ui - d x u 2 + d x 4> 2 

+ J J £ d x <p 2 dx - J p £ uid x (j) 2 dx - e j p £ u 2 d x <j) 2 dx. 



Using the L 1 uniform bounds for p £ and J £ , as well as the various Lipschitz 
bounds on (u±,u 2 ), it is clear that 



e J p £ u 2 d x <j) 2 dx + e j f £ {u\+eu 2 — v)dtu 2 dvdx 



< Ce. 



In the end, we have 

h 



i //.(«. 



+ eu 2 - v) 



d x ui + d x <pi 



dv dx 



+ / f £ (ui + eu 2 - v) d t ui + vd x ui - d x u 2 + d x <f> 2 



(2.34) 



dv dx 



+ J J £ \ -d x (/>i + d x (f> 2 ) dx 

- ^ J p £ ui(d x 4>i + ed x (j> 2 - d x 4>) dx 

- ^ J p £ u 2 d x ((f>i + e<t> 2 ) dx - j p £ E £ u 2 dx 

- j f £ \ui + eu 2 - v\ 2 d x ui + O(e), 

where 0(e) is a notation for all the terms which can be bounded by C £ , with 
C > independent of e. 
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Study of I2. 

We have: 



h = e j d t {d x 4>i + ed x (f)2){-d x 4> £ + d x 4>i + edxfa) dx 

-e j dtd x (/) e (d x (j>i + ed x 4> 2 ) dx. 
We get the easy bound (using \ab\ < ^(a 2 + b 2 )): 
£ J (d x (f>i + ed x (t)2){-d x 4> £ + d x <f>i + ed x (f> 2 ) dx 

<C e J {d t {d x 4>i + edxfo)) 2 dx + e j (-d x (j> e + d x (j>i + ed x <p 2 ) 2 dx 
< 0{e) + CU £ {t). 

Study of I3. 

We get: 

h = J d t (4>i+E<j)2)(-4>e + 4>i+£<p2)dx 

- J dtMfa + dx := /3 1 + II 
Let us start with l\, that we can rewrite as: 

J 3 = J (01 - 4>e)dt4>\ + £ j d t 4> 2 (-4>e + 4>i+ £<t>2) dx + e j d t <t>l(/>2 dx. 
If we differentiate with respect to time the Poisson equation, we obtain: 

-d t (f>e = -ed xx dt<f> £ ~ \ d tPe- 

Then, using the local conservation of charge (in shifted variables) that we 
recall below, 

--d t pe = — \d x p £ + -d x J e , 

we obtain that: 



^2 J d x p e (4>i + efo) dx 



+ \ j d x J £ (pidx + j d x J £ 4>2 dx 



- £ J d t dl x (j) e {<j)i+£<j)2). 

Of course, by integration by parts, we can rewrite the last three terms as 
follows: 

d x J e (pi dx + J d x J £ 4>2 dx + £ j dtdl x (j> e ((/>i + efa) 

= -\j JA^ dx - j J M 2 dx + £j d t d x Md^ + eOM 
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which get simplified using some terms coming from the computations of l\ 
and I 2 . 

Study of the remaining significant terms. 

We now have to study the following terms (coming from the computation 
of/i + / 2 +/ 3 ): 

Ki := - J p £ U2d x (4>i + e<j) 2 ) dx, K 2 := - J p £ E £ u 2 dx, 
K 3 := ~ J p £ ui(d x 4>i + ed x fc - d x <f> e ) dx, 
K4 := J 9xPe{<t>i + £<t>2) dx, 
K 5 := J (4>i - <p £ )d t <pi dx. 

We start with K\ and K 2 . Using the Poisson equation, we have: 

Ki = — J u 2 d x (<f>i + E(p 2 ) dx - e j <f> £ u 2 d x ((f)i + e<f> 2 ) dx 

+ e 2 J dl x 4> £ u 2 d x {4>i + e<j) 2 ) dx 

= J d x u 2 (j>i dx - e j 4> £ u 2 d x 4>i dx + e 2 J d 2 x 4> £ u 2 d x (4>i + e<f) 2 ) dx 

+ 0(e). 

We have (after integration by parts), 

e 2 J d 2 x 4> £ u 2 d x (4>\ + e(f) 2 ) dx <Ce. 
Similarly, we compute 

K 2 = j u 2 d x <j) £ dx + e J 4> £ u 2 d x <j) £ dx - e 2 j d 2 xx ^ £ u 2 d x ^ £ dx 

= - J d x u 2 <p £ dx - e J ^-d x u 2 dx + e 2 J d x u 2 dx. 

As a result, one gets 

K± + K 2 = y* (0i - <t> £ )d x u 2 dx + 0(e). 
Concerning K3, we have (once again using the Poisson equation): 

K 3 = ~ J ui(d x 4>i + ed x <j) 2 - d x (p £ ) dx 

- J 4> £ ui(d x 4>i + ed x (j) 2 - d x (p £ ) dx 

+ j edl x 4> £ ui{d x (j)i + ed x (f) 2 - d x <p £ ) dx 
:= Kl + K 2 + Kl 
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For i^g, we have: 



K\ = — j u\d x 4>i dx — J uid x (f>2 dx H — J u\d x (f> £ dx. 
Concerning K 2 , we write: 



u\d x u\ dx + / 4> £ d x (j) £ u\dx + e j § £ u\d x <\>2 dx 



and 



- / <P £ u ldxUl dx = - ^u ldxUl dx - j <p lUl d xUl dx. 

The contribution coming from j <j) £ d x (f) £ u\dx could be harmful (a priori 
it is of order 0(l/e)), but we can rely on the following identities in order to 
make the relative entropy appear: 



. , j) £ d x uidx 



j 4> £ d x <t) £ ui dx = - : 

= ~\ j ^ e ~ £( l )2 ) 2 ^ xUl dx ~ J ^e^ 1 ^ - " 1 dx + e J <p £ (j)2d x ui dx 

+ \ J^ 1 + £ ^ 2dxUl dx 
= ~ J (<Pe ~ 4>i ~ ^4>2) 2 d x ui dx + J (4>i - <p £ )<pid x ui dx + e j 4> £ ^2d x u\ dx 

~ \ J ^ xUx dx + e j ( < t > ifa d x u i + ^ £ 4>2 d xUi\ dx. 

We note that e J 4> £ 4>2d x ui dx is of order e, using the L 2 uniform bound on 
4> £ granted by the energy. 
Finally, for iff, we have 



d 2 xx 4> £ u\d x <\> £ dx + e j d 2 xx ^ £ u x {d x 4>\ + ed x ^ 2 ) dx. 



As before, the most significant term can be rewritten as follows: 



- j dl x 4> £ uid x (f) £ dx = e j ^^—d x ui dx 

= \ j £ ( d x4>e ~ d x 4>i - ed x ^ 2 ) 2 d x ui + e j d x <j) £ (d x 4>i + ed x 4> 2 )d x ui dx 
£(d x (fri + ed x 4>2) 2 d x tii dx. 



We have the bound: 



Jd x Ud^ + ed M d xUl dx 
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We treat K4 exactly as K3. 



Ki = ^ J p £ {d x 4>i + ed x (f) 2 ) dx 

= ~2 ji 9 ^! + ed x 4> 2 ) dx + ^- j 4> £ {d x 4>i + ed x (t) 2 ) dx 



J 



dl x (t> £ {d x (f)i + ed x 4> 2 ) dx 



:= K\ + K\ + K.\. 
Clearly, we have Kl = and 



4> e d x 4>i dx+ <p e d x 4>2 dx, 



and since J 4>id xxx (f)i dx = 0, we have: 



(f) £ )d xxx (f)i dx — el d xx (f> £ d x 4>2 dx 



By integration by parts, it is clear that the second term of K\ is of order 
y/i. 

Conclusion. 

Gathering all pieces together, we obtain that: 

f £ (ui +eu 2 - v)-{-d x u 1 +d x (pi)dvdx 



US) < n e (t) + 

+ / f £ (ui +EU 2 -V 



d t ui + u\d x ui - d x u 2 + d x (j) 2 



dv dx 



1 



+ - (ui- 4){)d x (j) e dx 



+ J (01 - 4>e) (d t (j>i + <t>xd x ui + uid x ui + dl xx 4>i + d x u 2 - d x 4> 2 )) dx 

- j(ui- 4>i)d x (j)2 - J (p\uid x ui dx - ~ J <\)\d x u\ 

- J f g \v - u x - eu2?d xUl dv dx-\ej - d xh - ed x ^d x 
) £ - 4>i - e§2) 2 dxU\ dx 



: u± dx 



ds. 



We impose the following cancellations (to kill all singular terms): 
' U\ - 01 = 0, 

(2.35) < d t ui + uid x U! - d x u 2 + d x <p2 = 0, 

k d t (pi + <t>id x u\ + u%d x ui + dl xx <j>i + d x u 2 - d x 4>2 = 0. 

These are consequences of the identity (12.221) . Using the Lipschitz bound on 
u\ we end up with: 

(2.36) U £ (t) < U £ (0) + !\cxH e {t) + C 2 ^) ds, 

Jo 
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which then yields the claimed Gronwall type bound f)2 .25[) . The proof is 
consequently complete. 

3. From the Vlasov-Poisson equation to the 
Zakharov-Kuznetsov equation 

We perform the analysis in 3D, but this can also be done in 2D, in an 
almost similar way. Let (ei, 62,63) be an orthonormal basis of IR 3 ; for sim- 
plicity we fix b = e\. 

We study the behaviour, as e — > 0, of the solutions to the following equa- 
tion (for (x, v) G T 3 x 



' edtfe - dxje +ev V x f £ + \E £ + ^5 ) • Y r f : - 0. 
E £ = -V x (j) £ , 



(3.1) 

e 2 A^ £ + e e ^ = / f £ dv, 
Jm.3 

k fe,\t=0 = fefi 

3.1. Preliminaries. This system possesses an energy which is conserved, 
at least formally: 

£e(t) '-=\ j fe\v\ 2 dvdx 

(3.2) 

+ 2 e / |V*&| 2 dx +~2 ( e - !) + !) 

Note that the third term of this energy has a L log L structure: 

i y (e £ ^ (e0 e - 1) + 1) dx = i y (e*** log(e* /e°) - e e ^ + e°) da;. 

We have the following global existence theorem, which can be adapted 
from the work of Bouchut [9j: 

Theorem 3.1. Let e > 0. Let f £j0 £ L 1 n L°°(T 3 x M 3 ) be a non-negative 
function such that: 

4(0) := -= I f £ ,o\v\ 2 dvdx 
(3.3) J 



+ 2 £ 



l -e j | V x £ ,o| 2 dx + ^ y (e £ ^-° (e</> £ , - 1) + 1) dz< +oo, 
where the initial electric potential (f) £ o is given by the elliptic equation: 
-£ 2 A^ e , + e e ^-° = / / £)0 dw. 

JK 3 



We a/so assume that: 



J f £fi dvdx = 1. 



TTien i/iere exists a non-negative global weak solution f e G Lf^^L 1 nL°°(T 3 x 
M 3 )) to ([33]) . suc/i i/iai: 

(3.4) Vt > t', £ E (t) < £ £ (i'), 
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and such that the following conservation laws for [p £ := J f e dv, J £ := J f e v dv) 
are satisfied: 

(3.5) d t pe - -d Xl p e + V x ■ J e = 0, 



(3.6) dt J £ - ~ e d Xl J £ + V x :(Jv® vf £ dv^j 



- -V^e^) + eV x ■ (V x V e V X V £ ) - £ -V x \V x V £ \ 2 . 

& — 

Let us now turn to the Cauchy problem for the ZK equation, that we 
recall now: 

(3.7) dt<f>i + <PidAi + d xl A Ul = 0. 

The only result in 3-D we are aware of for this equation concerns the case of 
the whole space M 3 (for results in 2D, we refer to |19[I34|): in |35j . Linares and 
Saut proved that ZK is locally well-posed in H S (M. 3 ), for s > 9/8, and more 
recently in [39], Ribaud and Vento showed that it is well-posed for s > 1. 
Their proofs are based on dispersive effects and can not be directly applied 
to the case of the torus T 3 . Using standard methods, we can nevertheless 
prove the easy theorem: 

Theorem 3.2. The ZK equation is locally well-posed in H S (T 3 ) for s > 5/2. 
This will be sufficient for our needs. 

3.2. Rigorous convergence result. Contrary to the ID case, we will not 
present the formal analysis allowing to guess that the limit equation is ZK. 
The principle is indeed identical, but the computations become quite lengthy. 
Let us refer to |29| for that point. 

We state directly the theorem asserting the convergence to ZK: 

Theorem 3.3. Let (/ e ,o)ee(o.i) ^ e a family of non-negative initial data sat- 
isfying the assumptions of Theorem Iff. II and such that there exists C > 
with: 

(3.8) S e (0) < C, Ve G (0,1). 

We denote by (/e) E e(0,i) a sequence of global weak solutions to (|3.1|) given 
by Theorem \3.1[ Let T~L e be the relative entropy defined by the functional: 
(3.9) 



n £ (t) :=\j 'fe 



I (1) (2) i 2 

4-1,, /£„M «,( 2 )| 2 _lU, /^W <rJ 2 )|2 

~T \V2 — y EUy — | + |t>3 — V £u 3 ~ £- u 3 I 



dv dx 



4>e l g( e £ (^) l^^+^+e 2 ^)^ _ e e4>s _|_ e e(<j>i+£(/>2+e 2 4>3)^j ^ 



1 

7 2 
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where 

(rh ,|W J 1 ) „W ,/ 2 ) J 2 ) ,/ 2 ) rh rh \ 

{<P\,U X ,U 2 ,U 3 ,U X ,U 2 ,U 3 ,(p2,(P3) 

G [C([0,T [,^+ 2 (T 3 ))] 2 x [C([0,T [,H s+1 (T 3 ))] 2 x [C([0, T [, H S (T 3 ))] 5 

(with s > 3/2 + 1, Tq > 0) satisfy the following system on [0,Tq[: 
(3.10) 

{ 2d t cpi + 2fad a .<t>i + d Xl (A + A ± ) Ul = 0, 

& = d* 1Xl </>i, 

4 2) = Cx 2 0i> 4 2) = d l^4>u 

(2) ^ 

d Xl u\ = d t ui + u\d x u\ + ^ ia . ia!1 0i, 
k <9zi</>3 = -d t (/)2 - d Xl {uf Vi) - «i 0*102 - u^d X2 4> 2 + U^'d X3 <j)2. 
Then there exist C±, C2 > 0, suc/i i/iai /or a// e G (0, 1), 

(3.11) Vi G [0,T [, H £ (i) < / (Ci-H e (s)ds + C 2 VF)ds. 

JO 

Assuming in addition that there exists C3 suc/i i/iai /or any e G (0, 1); 

(3.12) n £ (o) < c 3 ^, 

£/ien we obtain for any e G (0, 1): 

(3.13) Vt G [0,T [, ft e (f) < C 3 e Cl ^ + g2 eCl / l~ 1 ^- 

In this theorem, we need a large number of correctors in the relative 
entropy, precisely because of the nonlinear exponential term in the Pois- 
son equation (compare with the case of Theorem 12. 3| where this term is 
linearized). In order to get functions satisfying (|3.10p . we can proceed as 
follows: 

• The existence of <j)\ is ensured by Theorem 13.21 (actually the first 
equation of (|3,10p is slightly different from (|3.7p . but we can come 
down to (|3.7p by using some standard change of variables, see for 
instance |35]). 

• We observe that the six correctors ( \ ,62) have 
their value which is uniquely determined (contrary to the ID case). 
We accordingly set: 

( U W = ^ U W = -d^fa, ujp = d X2 4>i, 

02 = d XlX1 0i, 

(2) s 

^ d Xl u\ = d t ui + u\d x u\ + d XlXlXl <l>u ■ 

• Finally, <f>% is a high order corrector whose value is imposed by the 
last equation of (|3.10p . 
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Remark 3.1. Note that a global well-posedness result for ZK in T 3 would 
yield global in time estimates (that is Tq = +00 ) in this theorem, which would 
be significant for logarithmically growing times, as for Theorem \2.3l 

Remark 3.2. From (|3.13p . we deduce that: 

(3.14) J_ J ^l o g( e ^(0 e )/ e ^Wl+^2+e 2 03)) _ e e4>e + d% 

Following the terminology in the Boltzmann literature (see for instance the 
book |41| ). this roughly means that e £( ^ e "converges entropically" to e e (<t>i+ £ <f>2+ e <t>3) _ 

From the elementary inequality (|3.17p (which will be given later), we can 
also deduce the control: 

(3.15) ^ J (e^ - ebi^+^+^^y dx < C 3 e Clt ^ + C 2 1 yfl. 

We have as before the corollary: 

Corollary 3.1. Making the same assumptions as in the previous theorem, 
we obtain the weak convergences: 



(3.16) 



p e ^ £ _s.o 1 in LfM 1 weak-*, 

J £ -> e ^ (01,0,0) in LfM 1 weak-*. 



Up to some obvious modifications, the proof of Corollary 13. II is similar to 
that of Corollary 12.11 and therefore we omit it. 



3.3. Proof of Theorem 13.31 Relying on the fact the energy £ e (t) is a non- 
increasing function of time (and that T~L e is built as a modulation of the 
energy), we have for all t G [0, Tq[ and all e G (0, 1): 



Ue(t)<HM + f f 9t 



r / 1 l (1) , (2)|2 / (1) , (2)- 



+ 



+ 



+ e 



t 



dv dx ds 
dv dx ds 



d t 



d. 



e Jo 



f E + «4 2) i 2 - ^ (v^4 1} + «4 2) )) 

f £ (^\Vzu^ + evS^ 1 2 — V3 {\fiu^ + eu^)\ dv dx ds 

2 - V x Cp £ ■ (Va-01 + eV x fa) 
e £ <t>£ l g(1 /g^i+^a+e 2 ^)) _|_ e s(<t>i+£<t>2+s 2 4>3) 



-|V a .</>i + eV x <p 2 



dx ds 
dx ds 



:= H e (0) + / (II + I\ + If + I 2 + h) ds. 
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The general strategy in the proof will be to keep (without making approx- 
imation) all dangerous modulated terms of the form 



y/eu^ + eu^ -v 2 | [•••] dv dx 




and 

' (-e 6 * 8 + e e ^ +e * a+ea *'))[...] dx, 



where [...] contains some expression independent of e, as soon as a < 0. 
Then (|3.10|) is precisely designed so that all terms exactly vanish in the end. 

On the contrary, for a > 0, these terms will be of order e a (and thus 
small). For the first type of terms, this can be seen with the same argument 
as in the proof of theorem 12. 3| namely the uniform (in e) bounds on the L 1 
norm of p £ and J £ . For the second type of terms, one has to use the Poisson 
equation satisfied by <p £ and use the bound on the electric energy. Indeed, 
given some smooth function ^f, we can write: 



^dx = e a I e 2 A(j) e Vdx + e a I p £ ^dx. 



The second term is clearly of order e a , using the uniform L 1 bound on p £ . 
On the other hand, by integration by parts, we have for the first term 



e a J e 2 A(f) £ y dx = e a+2 J V4> £ ■ V* dx, 



which is of order e a+1 using the uniform bound obtained thanks to the 
conservation of energy: 



e J \V(f) £ \ 2 dx < S e (0) < C. 



Finally any term like e a f e e ^ 1+e< f' 2+e2 ^ [...] dx is clearly of order e a . 

We start by studying separately l\+ 1 2 + 1% and I3. All computations 
are justified using only the local conservation laws (13.51) and (|3.6|) . 



Study of /{ + If + 
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With similar computations as in the proof of Theorem 12,31 we obtain the 
identity 

l\ + l\ + If = 



(1) , (2) 
U\ + Eli) — V\ 



(2) 



dv dx 



\V^ U 3 } + £u 3 

- J ^PeU^ ■ V x 0i dx - J p £ u (2) • V x (j)i dx + J ^ J e • V x <pi dx 



\ 4 +^3^1 



7!"- 



y eu 2 + £ii 2 — ^ 2 



,( 2 ) 



-It. 



y/euy + eu\' - r 3/ 



,(2) 



(2) 



(if 



+ //. 



(2) 

x + eu\ - v\ 

(!) J_ =.„.( 2 ) 



+ £n (2) 



W + u^a^w - + a 



V 



-d xi uf> + 5x3 02 



X!<P2 



du dx 



+ j J £ ■ Vx02 dx — j p £ u^ ■ Vx02 dx — e J p £ u^ ■ V x (f>2 dx 

- J +euf ) -vxfd^ dx 

- j ^PeE e ■ u {1) dx- J p £ E £ ■ u (2) dx 
+ 0(e). 

In this equation, 0{e) is as usual a notation for all terms that can bounded 
by Ce, where C is a positive constant. We also denote here: 

U (D : = ( u U,^i u U,^ u W) and u( 2 ) := (u< 2 \ v%\ uf). 
We observe that we can write: 

- j -^PeE £ ■ dx - J ^PeU^ ' Vx01 dx - J p s U^> ■ V x (p2 dx 

= ~\ j Pz Ul ' ( V ^'01 + eV z02 - Vx0 £ ). 

Study of I2. 

We obtain: 

I i = £ f a(Vx0i+eVx02+e 2 Vx0 3 )-(-Vx0e+Vx0i+eVx02+e 2 Vx0 3 )dx 

- e y d t V x 4> £ ■ (Vx0i + eVx02 + £ 2 Vx03) dx. 

We have the easy bound: 
e / ^(Vx^i + eVx0 2 + £ 2 Vx0 3 ) • (-V x e + V x 0i + eVx0 2 + e 2 V x 3 ) dx 

< C e j (a t (Vx0i + eVx0 2 + £ 2 Vx0 3 )) 2 dx + e j {-d x <j) £ + d x <k + e^) 2 dx 
<0{e) + CU £ {t). 
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Study of I3. 

We can compute: 



s 



4 + 



The first term can be recast as follows: 
I l = \ J (-e e0E + e £( * 1+e<fe+£2<fe) )d^i dx + J (-e £ ^ + ^+ e ^+^ t M)d t <h dx 

+e f(-e £ ^ + e s <* l4 ** 24 * a * s >)fl^ 3 . 



For the second one, we use the Poisson equation (with a time derivative) 
and the conservation of charge in shifted variables, that we recall below: 

dt(e £ ^) = e 2 Ad t ^ + d t p £ , 

-d t pe = \d Xl p £ - -V x • J E , 
to obtain: 

+ \ j Je ^ X dX + j^ x ' Je ^ 2 dx + e j Je4>3 dx 



e j 9 t A0 e (^>i +e4> 2 + e 



2j. \ 
3j 



As in the proof of Theorem 12.31 every term in /| but the first one gets 
simplified with some other ones obtained in 1% and I 2 - 

Study of the remaining significant terms. 

There remain to study the following potentially harmful terms: 

Kf.= -J Pe uW ■ V x (<h. + e<fo) dx, K 2 :=- j p £ E £ ■ dx, 
K 3 := ~ J p e u W ■ (V x (pi + eV x 4>2 - V^ E ) dx, 
K 4 :=-^ J d Xl p £ ((f>i + e<j)2 + e 2 <j) 3 ) dx, 



\ [ {-e £ ^ +e £ ^ 1+£4,2+£2 ^)d t cp l dx. 
K 6 := I {-e £ ^ +e £ ^ 1+£ ' t>2+£2 ' t, * ) )dt<p2dx. 



K 5 : 

e 
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We start with K\ + using the Poisson equation: 

K l +K 2 = - j e £< M 2) • V x (<£i + e^2 - <f>e) dx 

+ e 2 J A0 £ uW . (V^i + eV x ^2 - V x e ) dx. 



= L x + L 2 . 



We have: 



Li = - y e e< M 2 ) • V^i dx-e J e £ ^n (2) • V x <k> 
For L2, we have 



-^2 + -^2- 



By integration by parts, we get: 



\Ll\ 



e 2 j V0 e • n (2) (A0i + eA0 2 ) dx 

<c[e 2 j \V<p £ \ 2 dx + Ce^j 
< C(e£ £ (t) + Ce) < Ce. 



Note here that we have used the Lipschitz bound on the second derivative 
of </>i and 02- 

For L 2 , we rely on the usual trick to write: 



= X -e 2 J \V x <f> e \ 2 V x -u<$ dx, 



from which we deduce that 



\L 2 \ < Ce. 



Using the Poisson equation for K3: 



K 3 = -^j e e ^uW . (V^i + eV x <h - V x £ ) dx 

+ J sAfeu^ ■ (V x ^i + eV x <f>2 ~ Vx<Pe) dx 
:=Kl + K 2 . 
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Concerning the first term, we write the decomposition: 

K\ = -- £ j e £ ^u^ ■ V x <Pi dx 

- J e^uW • V x (f> 2 dx + i J e £ ^u (1) • V x <f> £ dx 

:= Ji + J 2 + J 3 - 

We start by the study of J±: 
Ji = \f - e £ ^) uPd^fa dx 

~\ f e^'^u^d^ dx - i j eF*<{uW - u^ei) • V^i dx. 
Clearly, we have (for instance using a Taylor inequality): 
-i j e^-^u^d^cf)! dx = -^J u i l ) d Xl (f> l dx + j fau^ d^fa dx + 0{e). 

In the end, we will take = (pi, so the terms of this latest expression 
which are of order 0(l/e) and 0(1) are exactly equal to 0. 
For J2, it is sufficient to write: 

■h = ~ j e^u^d^fo dx- J e e ^(n (1) - u^ei) • V x fa dx. 
We decompose the last term in the following way 
■h = \f e £ ^4>id xl <p E 

+ iy e e +<{uP-<h)d xl <f> e 

+ iy e £ ^(u^ -u^ei)-V x( p £ . 

The first term could be dangerous, but will disappear using a term coming 
from a term of K4 (first term of K\ below). The second one will be equal to 

since we take = <pi. 

Gathering the pieces together, we finally obtain: 

X - j e e +°(uW - u^ei) ■ V x (4> £ -<Pi- e<p 2 ) = 

-4ft f ^(d^ + d xs u ( i ) )dx 



-±= j e £ ^ (u^d X2 (Pi+u^d x . A( Pi) dx 

- y/i J e £ ^ (u^d^ + u^d^) dx. 



Let us now turn to the treatment of K\ : 
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This term is treated exactly as L2, but we have to be careful since it is 
singular in e. Here, rather than bounding by the energy, we shall rely on a 
bound by the modulated energy. We have 



= -e/2 J \V x cP £ \ 2 V x -u^ dx 
- e/2 j \V x <f> £ - V^i - eV x (f>2 - e 2 V x <p 3 \ 2 V x ■ dx 
+ e/2 J \V x <f>i +eV x <f> 2 + e 2 V x <f> 3 \ 2 V x -u^dx 

- e J V x <j) £ ■ (V x 0i + eV x fa + £ 2 Vxfe)V, • u {l) dx. 

We shall bound the very last term using the Cauchy-Schwarz inequality 
and the various Lipschitz bounds: 

e j V x <f> e ■ (V x <fc + eV x <h + e 2 V x (j) Z )V x ■ «W dx 
< CVlJe J \V x <f> e \ 2 dx < CVe~£ £ (t) < C^I. 



Therefore we obtain: 



\k\\< cn £ (t) + c^i. 



Finally, concerning K4 we get: 

K4 = ^ J Pe{d Xl 4>\ + ed Xl 4> 2 + e 2 d Xl <j) 3 ) dx 

= ~^f z £<t>e (d xl (i>i+ed Xl (t)2 + e 2 d Xl <j)z)dx 

- j A(f> e (d xi <f>i +ed Xl 4>2 + e 2 d xi <j) 3 )dx 
:= K\ + K\. 
The first term can be decomposed as: 

Kl = ~J d Xl 4) £ e e ^4>i dx 
+-^f e e ^d xx <\>2 dx 
+ J e E ^d xl cf) 3 dx. 
The second one can be recast as: 

K\ = - j (/> £ Ad xl <f>i dx + 
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We shall focus our attention on the important term: 
L :=- J <p e Ad Xl (pi dx 

= ~\j e^Ac^i dx + ~ J [e £ ^ - s^Ad^ dx. 

We have the following technical result, allowing to consider the second 
term above as an error term: 

Lemma 3.1. There exists C > such that for any e G (0, 1): 
1 



J (e £ ^ - e<j) E ) Ad xi 4>i dx 



Proof of Lemma \3.1\ The naive idea would be to use the Taylor expansion 

e x ~ 1 + x + ^x 2 . 

but we cannot say that e||</> e ||oo 1- (Even worse, we do not have any 
L 2 control on cf) £ .) Instead, we will only rely on the bounds given by the 
conservation of energy. The classical inequality (valid for x,y > 0) will be 
very useful: 

(3.17) {y/x - y/y) <x\og(x/y) - x + y. 

We shall write the decomposition 

\ j (e £ ^ -e</> e )Ad xl <l> 1 dx = 

i J (e e ^ - 2e^)Ad Xl 4>! dx 

+ ~ J(e^ -^e^Ad^dx. 
We first recast the fist term as follows: 

i J (e £ ^ - 2e^)Ad xl fa dx = J (e £ ^ - 2e^ + l)Ad Xl <f>i dx 

= \J ieh £ ^ - ifAd^fa dx. 

Using (|3.17|) and the Lipschitz bound on the second order derivative of 
<pi, we obtain: 

I J ( e §#* _ ifAd^fa dx<Ce j ^{e £ ^{e4> £ - 1) - 1) dx 

< eS £ (t) 

< e£ £ (0) 

< Ce. 

For the second term, we have by integration by parts: 
I J ( e |^e _ I^Ac^i dx = - J VxMe^ ~ 1) • VA^i dx. 
Then we write, using \ab\ < ^(Xa 2 + jb 2 ), with a a parameter to be fixed: 
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1) • Vd Xl 4>i dx 



KCe^e J |V^ £ | 2 + Ce~ a+2 J ^(e £ ^ (e^ - 1) + 1) dx. 

To optimize, it is clear that we should take a = 3/2. Thus, relying on the 
uniform bound given by the energy, we get: 



e2 



2Ad xl (pi dx 



< Cy/e. 



The proof is therefore complete. 



Conclusion. 

Finally, we impose all cancellations for the terms of the terms 

/ (1) , (2) 



£ 



fe 



v 2 | [•••] dvdx 

v 3y 



for a 



-1, -1/2,0 and 



for j3 = —3/2, —1, —1/2,0, which precisely means that we impose: 
(3.18) 

- d Xl vfp + d Xl 4>i = 0, 

- 4 + d X2 (f)i = 0, 
u)p + d xs (f>i = 0, 



u 



(2) 



d xl U 2 



^2 d Xl Ur^ 



0, 

,(!) 



(9 t u^ + u\ X) d Xl u\ ±J - d Xl u { 2 l) + 9^02 
0, 
0. 



d xi v£p + 3^02 



^2 ^ x 3 ^3 

,(2) 



0, 



□ 



d t (j)i + V x ■ vS z > + fad^fa + A9 Xl ^i - ^2 = 0, 
w^^^i + u^d X3 (pi = 0, 

k dt4>2 + 4>\d Xl 4>2 + uf'dx^i + v§'d X2 4>x + uf'dx^i - d Xl 4> 3 = 0. 

All these indentities are exactly obtained as a consequence of (|3.10p . 
We end up with: 

(3.19) n s (t) < H £ (0) + / {CxHS) + C 2 Ji) ds, 

Jo 

and arguing as in the proof of Theorem 12.31 this is over. 
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4. Appendix A: From the Vlasov-Poisson equation to the 
Kadomstev-Petviashvili II equation 

Assuming a slow variation in the #2 direction, one may end up with the fol- 
lowing anisotropic long wave scaling for the Vlasov-Poisson system for ions. 
For simplicity, we restrict here to the linearized Maxwell-Boltzmann law, but 
the same study can be performed for the full equations. The variables are 
t > 0,x £ T 2 ,v G R 2 : 

' ed t f £ - d Xl U + e Vl d Xl f £ + s 3 / 2 v 2 d X2 f £ + E £ - V v f £ = 0, 
E = (-d xi <f> £ , -y/ed X2 (f) £ ), 



(4.1) 



2f)2 



e z d; 



3a2 



e A d: 



+ ec 



fedv - 1, 



^ f\t=Q = fo- 

The scaled energy of this system is the following functional: 



(4.2) 



£e{t) := \ j fe\v\ 2 dvdx 

+ \ £ J \d Xl (f> £ \ 2 dx + ^e 2 J \d X2 </) E \ 2 dx + ^ J 



4> 2 dx. 



We have the existence of global weak solutions, sharing the same properties 
of those given in Theorem 12.11 
For the KP-II equation, that is 

(4.3) d Xl {d t 4>i + <t>id xi 4>i + dl 1XlX1 4>i) + d 2 x ^ x = 0, 

our reference is an article by Bourgain |12| . in which is proved the following 
theorem: 

Theorem 4.1. The ZK equation is globally well-posed in H S (T 2 ) for s > 0. 

Once again, we will only use this theorem for large values of s. 
As for the other cases, we obtain the rigorous convergence to KP-II, which 
is summarized in the following theorem: 

Theorem 4.2. Let (/e,o)ee(0,i) ^ e a, family of non-negative initial data such 
that 



(4.4) 



l|/e,Q||£ 1 nL°° < +°°) J fefldvdx = 1, 



and such that there exists C > with: 

(4.5) S e (0) < C, Ve G (0,1). 

We denote by (/ e ) £ g(o,i) a, sequence of global weak solutions to (|4.ip . Let 
T~L £ be the relative entropy defined by the functional: 
(4.6) 



n £ (t) :-- 



1 



fe 



I (!) 



eu 



(2) ,2 



+ \V2 - V^U 



(1) 

2 



3/2 u (2) |2 



dv dx 



+ \ £ j 0e ~ 3bi<^i - edxifal 2 dx + ^e 2 J \d X2 4> £ - d X2 4>i - ed, 

">e ~ 4>l ~ £ 4>2) 2 dx, 



X2<P2\ 



1 

+ 2 
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where u 2 , fa, fa) € [C([0, +oo[, H S+2 (T 2 ))] 4 (with s > 2) satisfy the 
following system: 

( d Xl (2d t fa + 3fad Xl fa + d 3 XlXlXl fa) + d 2 X2X2 fa = 0, 



(4.7) 



,( 2 ) 



'I, 



5x 2 



k d Xl (u\ ' - fa) = d t fa + fad Xl fa. 
Then there exist C\,C 2 > 0, such that for all e G (0, 1), 

(4.8) H £ (t) < He(P) + f {C 1 H £ {s)ds + C 2 ^)ds. 

Jo 

Assuming in addition that there exists C3 > such that for any e G (0, 1), 

(4.9) H £ (0) < C^e. 
Then we obtain for all e € (0, 1): 



(4.10) 



Vt G [0, +00 [, H £ (t) < C 3 e Cl Ve + C 2 - 



1 



Ci 



as well as the weak convergences: 

{ p E ^ e ^o 1 in LfM 1 weak-*, 

(4.11) r ' 

[ J £ ^e->o (01,0) in LfM weak-*. 

The proof of Theorem 14.21 follows from computations in the same spirit as 
the previous ones, and therefore we leave it to the reader. 

5. Appendix B: A KdV limit in the whole space R 

All results stated in this paper are restricted to PDEs set in the torus 
for the space variable. The reason is that in all cases, in the end, the first 
moment p £ (charge density) has to weakly converge to the constant 1. This 
function is obviously not integrable in the whole space, and the assumptions 
needed for our results to hold are actually not consistent in the whole space 
case. 

It is nevertheless possible to slightly adapt the KdV limit of Section [2] to 
handle thatcase. Keeping the same notations, we shall rely on the fact that 
the electric potential can be defined up to a constant, and use the following 
version of the energy 
(5.1) 



m) ■■= ~ / 



1 



f e \v\ dvdx + -e / |Va;0 e | dx + — 



1 



1 



i E ) dx dx, 



This means that fa £ - + L 2 (R) (instead of the more usual fa G L 2 (R)). 

The Ansatz in that case for the formal computations now corresponds to 
the following one: 

' Pe = ep 1 + 0(e 2 ), 
1 



(5.2) 



, £ - - + fa+efa + 0(e 2 ), 
l £ = ui + eu 2 + 0(e 2 ). 
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In particular this means that we expect that p £ weakly converges to (which 
is of course integrable on R). The formal computations then remain the same. 
In the end, we may obtain the same result as Theorem 12.31 except that we 
consider the following relative entropy instead of f|2.2 1 j) : 
(5.3) 

% E (t) := - / f £ \v-ui-eu 2 \ 2 dvdx + -e \d x <j) e - d x <f>i - edxfo] 2 dx 
* JRxR ^ Jr 

+ o/ ( 4>e <j>l - £02 ) dx. 
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